Stamping and wrinkling of elastic plates 



Jeremy Hure, Benoit Roman, and Jose Bico 
PMMH, CNRS UMR 7636, UPMC & Univ. Pans Diderot, 
ESPCI-ParisTech, 10 rue Vauquelvn, 75231 Paris Cedex 05, France. 

We study the peculiar wrinkling pattern of an elastic plate stamped into a spherical mold. We 
show that the wavelength of the wrinkles decreases with their amplitude, but reaches a maximum 
when the amplitude is of the order of the thickness of the plate. The force required for compressing 
the wrinkled plate presents a maximum independent of the thickness. A model is derived and verified 
experimentally for a simple one-dimensional case. This model is extended to the initial situation 
through an effective Young modulus representing the mechanical behavior of wrinkled state. The 
theoretical predictions are shown to be in good agreement with the experiments. This approach 
provides a complement to the "tension field theory" developed for wrinkles with unconstrained 
amplitude. 



Wrinkling patterns are observed when thin plates are 
put under compression, spanning scales from geologi- 
cal patterns [T], skin wrinkles resulting from aging pro- 
cesses or scars [3] to cells locomotion generating strains 
on substrates [3]. They have important applications in 
micro-engineering such as the formation of controlled 
patterns [5], or the estimation of mechanical properties 
from the number and the extent of the wrinkles [5HS]- 
While most studies have focused on near threshold pat- 
terns, recent contributions have pushed further the de- 
scription of finely wrinkled plates, i.e., well above the 
initial buckling threshold [TUl E] • These studies consider 
plates submitted to strong in-plane tension on the bound- 
aries, which prevents stress focusing commonly observed 
in crumpled paper [12] . In these descriptions, wrinkles 
are assumed to totally relax compressive stresses, as in 
traditional tension field theory [T31 E] • In addition, both 
amplitude and wavelength of the wrinkles vanish with the 
thickness of the plate. We propose to study a conceptu- 
ally different wrinkling regime where the boundaries are 
free from tension but the amplitude of the wrinkles is 
highly constrained. We focus on the simplest example, 
an elastic plate compressed in a spherical mold with a 
confinement defined by a gap 8 (Fig. [I]) . This stamping 
configuration is common in industrial processes where 
metal plates are plastically embossed, the mismatch in 
Gaussian curvature generally leading to regular wrinkles 
[T5HT7] . In the elastic case, crumpling singularities first 
appear as the mold is progressively closed down (Fig. [TJ;) 
and evolve into a pattern of apparently smooth radial 
wrinkles for high confinement (Fig. [Iji-e). We show that 
constraining the amplitude does not lead to the collapse 
of compressive stress. Instead, the wrinkling pattern 
derives from a nontrivial balance between compression 
and bending stresses, with surprising consequences : the 
wavelength of the wrinkles does not vanish and the con- 
straining force reaches a maximum independent of the 
thickness of the plate. 

One- dimensional problem. We start by considering the 
simpler problem of a plate of length L, thickness h and 
unit width. In-plane displacements are imposed at both 
ends u x (±L/2) = ±A/2 (Fig. [2). In addition, the out-of- 
plane displacement of the plate is constrained to a max- 
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FIG. 1: Top: Experimental setup. A circular plate of radius 
R is compressed between two rigid transparent hemispherical 
dies of radius p. The gap S between the spheres is imposed. 
Bottom: Experimental observations of wrinkles formed when 
a plate is compressed between two hemispherical dies. S is 
decreased from (b) to (d) while R is maintained constant, (d) 
and (e) are characterized by the same S, but different values 
for R (p = 58mm, h — 58 /im and R = 12.5/25 mm). 



imal value S. The amplitude A of the median plane of 
the plate thus corresponds to A = S— h. For large values 
of S a single wrinkle forms as the axial force exceeds Eu- 
ler's critical load [TB] ■ Within small slope approximation, 
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FIG. 2: Experimental setup: a plate of length L, thickness h 
and unit width is subjected to symmetrical displacements at 
the ends u x (±L/2) = ±A/2 and a maximal deflection 5. We 
consider that n similar wrinkles are formed (A „ e = A/n) of 
wavelength A = L/n and amplitude A, with A = S — h. 



the out-of-plane deflection of a single wrinkle (Fig. [2| , is 
taken as w(x) — A/2[l + cos (27rx/A)], where A is the 
wavelength. The associated bending energy thus reads: 
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with B = Eh 3 / [12(1 — v 2 )] the bending modulus of the 
plate (E and v the material Young's modulus and Pois- 
son's ratio). We are here interested in the non-classical 
limit 8 — > h where in-plane stress a and strain e can- 
not be neglected. According to the in-plane equilibrium 
equation d x a = 0, a and e are constant along the plate, 

leading to e = f^Jdu/dx + (l/2)(dw/dx) 2 } dx/X. The 
corresponding stretching energy can be written as: 
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where A one corresponds to the in-plane displacement for 
a wrinkle and S — Eh/ (1 — v 2 ) is the stretching modulus 
of the plate [2"5] , 

We consider a series of n successive sinusoidal wrin- 
kles along the length L of the plate. Such model does 
not reproduce the continuous evolution of the confined 
plate, as flat parts are actually observed [TH1 [5D] • Nev- 
ertheless, it describes exactly particular shapes taken by 
the plate all along the stamping process. To simplify 
the description, we thus propose to join these particu- 
lar states continuously by considering non-integer values 
of n. The global energy corresponding to n successive 
wrinkles is obtained by using the conditions nX = L and 
nA one = A in equations ([!]) and (|2|): 
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Minimizing the total energy for an imposed amplitude 
with respect to the number of wrinkles leads to: 
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FIG. 3: Maximal number of wrinkles for polypropylene films 
of thicknesses ranging from 15 fim to 250 /im. Open symbols, 
ID experiments: strips of length 250 mm, width 25 mm with 
A/L — 0.2% and 0.3%. Filled symbols, 2D experiments: discs 
of radius R ranging from 10 mm to 70 mm, radius of spheres 
p = 58, 600 mm and n ma x the number of wrinkles at the edge 
of the discs. For a comparison with ID experiments, we take 
L = 2tvR and A/L = R 2 /8p 2 [2T]. The solid line corresponds 
to the best fit of the data: n max = 0.20 (L/K)^J 'A/L. 



Counter-intuitively the number of wrinkles tends towards 
a finite value n rnax ~ 0.39(L//i)a/ A/L as the ampli- 
tude A vanishes [55] . In order to validate this prediction, 
polypropylene films (E = 2200 MPa and v = 0.4, Innovia 
Films) of thicknesses h ranging from 15 /im to 250 /jm are 
cut in bands of length L = 250 mm and 25 mm of width. 
The strip is clamped on a rigid plate with an imposed 
displacement A. The resulting blister is confined as illus- 
trated in Fig. [2] Note here that to avoid plastic events, all 
the experiments are such that A/L < ey, where ty ~ 2% 
is the yield point of the material. We observe that for 
large confinement the number of wrinkles saturates to a 
maximum value in good agreement with the theoretical 
power law (Fig. [3]). Nevertheless the value of the pref- 
actor is significantly lower than predicted (0.2 instead of 
0.39). We interpret this discrepancy as a consequence 
of the discrete transition between successive modes |19j . 
Indeed the buckling of n wrinkles can result into n + 1 to 
Zn wrinkles (if all wrinkles split simultaneously), which 
leads to an ambiguity in the actual mode number. 

In addition to the number of wrinkles, the load (per 
unit width) required for compressing the plate can finally 
be derived by differentiating the elastic energy: 
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where a = A/h is the relative amplitude. Within the 
limit A ^ h, both forces are proportional to the bending 
stiffness Eh 3 (the plate can then be considered as inex- 
tensible) and increase with the confinement. However the 
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FIG. 4: Evolution of the vertical force (per unit width) f z re- 
quired to confine the plate as a function of the dimensionless 
amplitude A/h, defined as f z (A/h = 0.94) = f z , ma x- The 
solid line corresponds to eq. |5b| the dashed line to the in- 
extensible approximation. Experiments with strips of length 
250 mm, width 25 mm with A/L = 0.07% (circles) and 0.2% 
(squares). Inset: Raw force data for h = 90 /im, solid line 
corresponds to A/L — 0, squares to A/L — 0.2%. 



vertical load vanishes for high compression (a <C 1) af- 
ter reaching a maximum f z , ma x — .2AE 'L(A / L) 2 (with 
v = 0.4) for A/h ~ 0.94, while the axial force (per unit 
width) tends towards (A/L) S. When A — 0, we recover 
the simple case of a flat plate under lateral compres- 
sion where buckling is inhibited, and f z vanishes. The 
experimental measurement of the vertical force is del- 
icate since any slight misalignment generates an addi- 
tional force during the compression. Capturing the limit 
A — > also requires an accurate positioning of the rigid 
dies. To limit artifacts due to misalignment, the force 
due to the wrinkles f z is obtained by substracting the 
force given by the tensile machine for A/L = (corre- 
sponding to a flat plate) to the one for A/L ^ 0. As 
shown in the inset of Fig. |4j this difference clearly indi- 
cates the presence of a maximum, as predicted by our 
description. The precise localization of A = is still dif- 
ficult and is inferred from the position corresponding to 
the maximum force, f z (A/h = 0.94) = f z , max (Fig. B. 
We observe that for a fixed value of A, the axial force in 



eq. (5a) is proportional to the applied strain A/L, which 



brings us to define an effective Young's modulus: 
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This effective modulus increases from zero to E as the 
plate is progressively confined. 

Two dimensional patterns. We consider elastic discs 
(same material) with radius R (from 10 mm to 70 mm) 
and thickness (from 30 /mi to 250 /mi) embossed between 
hemispherical stamps of radius p (58 mm and 600 mm), 
separated by a distance 5 (Fig. [Tj) . The compaction force 
is measured as the difference of the raw force with the 
force obtained without the plate (Fig. [5]). The orthora- 
dial compression responsible for the formation of wrinkles 



is a consequence of Gauss Theorema Egregium: wrapping 
a sphere with a planar sheet implies stretching or com- 
pression in addition to bending [25]. In terms of scaling, 
the perfect contact between the plate and the stamp in- 
volves a typical strain e ~ (R/p) 2 and a stretching en- 
ergy EhR 6 /p 4 [2TJ [23]. The ratio of the stretching en- 
ergy to the typical bending energy Eh 3 R 2 /p 2 thus yields 
the dimensionless parameter R/^fph. Therefore form- 
ing wrinkles is not expected to release energy in the case 
R 2 <C ph, i.e. when the bending is dominant over stretch- 
ing energy. In the experiment, if R is small enough, the 
plate indeed remains unwrinkled when the mold is closed. 

We focus on the opposite limit R 2 3> ph, especially for 
8 ~ h, where hierarchical wrinkles appear beyond a cer- 
tain distance from the center of the plate. The minute 
amplitude of the wrinkles lead us to assume that the 
wrinkles only involve smooth deformations as a pertur- 
bation from the state S = h, where the plate lays un- 
wrinkled on the spherical cap. In other words, we con- 
jecture [57] that the observed transition in fig. [TJ;,d cor- 
responds to stress defocusing similar to the transition de- 
scribed in [24]. Within this assumption, the scaling law 
for the maximal number of wrinkles derived in the ID 
situation can be adapted to the circular plate. Since the 
effective orthoradial strain is given by (r/p) 2 , the number 
of wrinkles observed at the edge of the plate for high con- 
finement is expected to follow n max ~ R 2 / ph. To com- 
pare quantitatively ID and 2D experimental results, we 
consider that L = 2wR and A/L — R 2 /8p 2 , which corre- 
sponds to the strain at the edge of a plate completely in 
contact with the sphere. Fig. [3] shows a good agreement 
between the prediction for the wavelength and the ex- 
periments (with the same prefactor as in ID). Note that 
using eq. [4] near the middle of the plate where egg = 
would give n = 0, which is not the case. The detailed 
description of the wrinkling cascade can not be captured 
by a simple ID model and the transition between flat and 
wrinkled parts of a plate has previously been shown to 
be problematic |10j . However our scaling would predict 
a number of wrinkles at the middle of the plate n max /A. 
If the transition were set by period doubling this number 
would correspond to three generation of wrinkles, which 
is qualitatively observed (Fig. [Tj) . 

In order to describe the extension of the smooth re- 
gion and the stamping force, we consider a model plate 
with an effective Young modulus E e f f in the orthoradial 
direction and E in the radial direction. As a first approx- 
imation we assume that wrinkles appear for agg < 0. We 
thus consider E e f f — E for agg (r) > when no wrinkles 
are formed and E eff = E/[l + § (A/h) 2 } for cr ee (r) < 0. 
The constitutive law for this effective anisotropic mate- 
rial can be derived from the initial Hooke's law: 
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where we take for simplicity the Poisson ratio of the ef- 
fective plate equal to zero. The corresponding strains are 
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given by e rr = du/dr + (l/2)(dw/dr) 2 and egg = u/r, 
with u and w are the radial displacement and the de- 
flection of the plate, respectively. In addition we as- 
sume the effective stress field in the finely wrinkled re- 
gion is axisymmctric a r g = 10J. Within the limit 
of high compression (5 ~ h), w is expected to follow 
w(r) = —r 2 /2p. The different strains and the location 
a of the transition between the smooth and the wrin- 
kled regions are finally derived by solving the equilibrium 
equation da rr /dr + (<j rr — ogg) jr = with the boundary 
conditions it(0) = 0, a rr (R) = 0, cr rr (a~) = a rr {a + ) and 
agg(a) — 0. The analytic solution provides a description 
of the confined plate from finite to zero amplitude. Al- 
though the radial stress is always tensile, the orthoradial 
stress is only tensile in the central region of the sheet 
and progressively becomes compressive towards the pe- 
riphery where wrinkles are observed (Fig.[T]). In the case 
of high confinement, the radius of the unwrinkled zone 
is a = i?/vo, in qualitative agreement with our experi- 
mental data. 
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FIG. 5: Evolution of the vertical force F z i required to con- 
fine the plate as a function of the dimensionless amplitude 
A/h, defined as F z \{A/h = 1.14) = F z i )max . Experiments 
with discs of radius 25mm, thicknesses h = 90 /im (squares), 
150 /im (circles) and 250 /im (triangles), with p = 58mm. 
The solid line corresponds to the theoretical prediction. In- 
set: Raw force data, the solid thin line corresponds stamping 
without plate, circles to h — 150 /im. 



with a spherical shell with the same material properties 
as the compressed sheet [28] . 

Nevertheless the evolution of the stamping force with 
the dimensionless amplitude exhibits a fair agreement 
with the prediction (Fig. [5]). In addition, we verified that 
the maximum of the force follows F z \_ max = [(1.85 ± 
0.05)l(T 3 ]£ , i? 6 / / 9 4 - The prefactor given by our model is 
4.10 -3 . The discrepancy is however consistent with the 
difference between the predicted and observed maximal 
number of wrinkles in ID: considering that the number 
of wrinkles is half the one given by eq. [4] would indeed 
decrease the estimate of the maximal vertical force f z by 
a factor close to 2. 

To conclude, we have studied experimentally and the- 
oretically the packing of wrinkled thin plates constrained 
in amplitude. We have broadened the one-dimensionnal 
analysis |19j for vanishing amplitudes with two novel 
characteristics: the number of wrinkles tends towards a 
finite value and the stamping force exhibits a maximum 
value. We have finally shown how these results can be 
relevant in a 2D geometry by using an effective plate that 
takes into account the presence of wrinkles. The effective 
compressive Young's modulus of the wrinkled zone which 
depends on the imposed wrinkle amplitude is analogous 
to the complete collapse of compressive stress in plates 
with unconstrained amplitude but under large tension. 
Our approach can be applied to other stamping geome- 
tries, such as negative curvature shapes where we expect 
different wrinkling pattern. 
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In addition to the description of the compressed plate, 
the combination of both ID and 2D approaches also pro- 
vides a simple estimate of the force required for stamp- 
ing the plate, which is relevant for practical applica- 
tions. Indeed reducing the amplitude of the wrinkles rc- 
quiresa force F z \ that was derived for the ID situation 
(Eq. 5b |. By replacing L by R and A/L by (R 2 /p 2 ), we 
thus expect F z i cx ER 6 /p 4 . More precisely, we compute 

this force as F z i = J f z 2irrdr (with A/L — egg and 
L = 2nr). Moreover, curving the plate induces an addi- 
tional pressure P = h(a rr + <Jgg)/ p and thus a curving 
force F Z 2 ~ EhR 4 /p 3 . However, in the relevant limit 
R 2 S> ph, we expect F Z 2 <C F z \. Measuring the actual 
force is more delicate than in the ID experiment since ob- 
taining a reference force would require a sphere covered 
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